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Abstract 


One of the earliest approaches in gain-scheduling control is the gridding based approach, in which a set of local linear 
time-invariant models are obtained at various gridded points corresponding to the varying parameters within the flight 
envelop. In order to ensure smooth and effective Linear Parameter- Varying control, aligning all the flexible modes 
within each local model and maintaining small number of representative local models over the gridded parameter 
space are crucial. In addition, since the flexible structural models tend to have large dimensions, a tractable model 
reduction process is necessary. In this paper, the notion of o-shifted #4- and #-norm are introduced and used as 
a metric to measure the model mismatch. A new modal alignment algorithm is developed which utilizes the defined 
metric for aligning all the local models over the entire gridded parameter space. Furthermore, an Adaptive Grid 
Step Size Determination algorithm is developed to minimize the number of local models required to represent the 
gridded parameter space. For model reduction, we propose to utilize the concept of Composite Modal Cost Analysis, 
through which the collective contribution of each flexible mode is computed and ranked. Therefore, a reduced-order 
model is constructed by retaining only those modes with significant contribution. The NASA Generic Transport 
Model operating at various flight speeds is studied for verification purpose, and the analysis and simulation results 
demonstrate the effectiveness of the proposed modeling approach. 
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1. Introduction 


At the advent of advanced composite materials technology, lightweight aircraft design concept has attracted con- 
siderable attentions in recent years in an effort to improve the aerodynamic efficiency. However, as the flexibility of 
a structure increases, aeroelastic interactions with aircraft aerodynamic forces and moments can have adverse impact 
on aircraft’s stability and performance. Therefore, active control of aeroelastic aircraft is becoming increasingly im- 
portant, especially for aircraft with highly flexible wings. The objectives of flight control design are to guarantee the 
closed-loop stability and to improve the handling qualities over a wide range of flight conditions. Generally, these 
objectives cannot be achieved using just one single flight controller operating over the entire flight envelop, due to the 
notable variations on rigid body aerodynamics and aeroelasticity within the flight profile. 

One effective remedy for this challenging problem is to directly synthesize Linear Parameter- Varying (LPV) con- 
troller [1] for the entire flight envelop [5]. This approach has been demonstrated to be very promising for 
various aerospace applications [6] [7| [8] {9}. Although many theoretical results were developed and applied to the spe- 
cial class of nonlinear and/or parameter-varying systems, the essential requirement for these theories to work is to 
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have a representative LPV model that captures the varying nature of the system dynamics. In general, there are two 
ways to attain the needed LPV model. If a mathematical model is difficult to attain [10], a global{11] or local (12) 
system identification technique can be applied to identify the LPV model. In the global system identification ap- 
proach, a single experiment has to be performed such that all scheduling parameters are excited simultaneously and 
persistently [13], which is not always possible in practice. On the other hand, local system identification approach 
requires performing multiple experiments with fixed scheduling parameters at various operating points, leading 
to a set of local linear time-invariant (LTI) models that need to be interpolated to obtain the LPV model. However, the 
main drawback associated with local methods is the lack of information for the rate of change of scheduling parameter 
since these LTT models are obtained by freezing scheduling parameters at predefined points. 

In the context of attaining aeroelastic models suitable for robust and LPV control, a number of research has been 
conducted recently, see for instance [15] [16] [T7] [18]. Some of these studies assume the availability of the analytical 
models{19] [21], which is not practical since most of the large scale aerospace structures are obtained via finite 
element analysis method, hence they are discrete in nature. In Varga et al. and Puyou and Losser [23], high- 
fidelity models were considered that included both rigid and flexible dynamics of a civilian aircraft. However, direct 
interpolation is not possible due to the inconsistency of their state-space representations at different operating points. 
Similar results were reported in Ferreres and Roos [24], where empirical criteria and physical knowledge of 
the system were required to obtain the LPV model. Recently, the problem of modal matching and model reduction 
were considered in Theis et al.[15]. A set of criteria were introduced for modal alignment, based on the modal 
frequency, damping ratio, and the ”B” and ”C” matrices associated with each mode. For model reduction, truncation, 
residualization, and balanced coordinate transformation approaches were utilized to reduce system order sequentially, 
and the v-gap metric was used to measure the distance between the reduced-order models (ROMs) and full-order 
models (FOMs). These approaches were applied to study the Body Freedom Flutter (BFF) air vehicle. Similar 
approach was used in for BFF aircraft structure. In , a procedure was developed to derive reduced-order 
LPV model from a given set of large-scale LTI models at grid points. #-norm was used to measure the distance 
between FOMs and ROMs. However, the approach was formulated in the Linear Fractional Representation (LFR) 
framework. The tensor product model transformation technique was proposed in Baranyi and Takarics [18], 
in which the LPV state-space models were transformed into parameter varying convex combination of LTI models 
through high-order singular value decomposition process. As a result, only those LTI models representing polytopic 
vertices were needed for subsequent control design. Although modeling and control of aeroelastic aircraft in LPV 
framework have been studied extensively in the literature, the topic still remains as a challenging research subject and 
needs to be explored. 

The main contributions of this paper are as follows. First, we reexamine the linear interpolation approach. It is 
found that when interpolating two aeroelastic LTI models operating at two different flight conditions, the resultant 
interpolated model can be erroneous, especially when performed in modal coordinates. This is because the sequential 
order of the flexible modes at two different operating conditions might be different, hence resulting in a model mis- 
match when interpolating the two models. Therefore, a special care must be taken to ensure that the flexible modes 
are consistently ordered throughout operating conditions before performing interpolation between any two LTT mod- 
els. For this purpose, a novel Modal Alignment Algorithm (MAA) is developed to align all aeroelastic modes for 
a given set of LTI models sampled over the flight envelop, so that all LTI models have consistent modal sequence. 
Second, the number of LTI models needed to approximate the aerodynamic variations within the flight envelop can be 
minimized. With the current advanced computational power and software technologies, it is not difficult to generate a 
large number of LTI models with very fine grid points that cover the entire flight envelop. While this large number of 
LTI models helps improve modeling accuracy, it unnecessarily burdens the flight control design and implementation 
efforts, since a flight controller must be designed for each local LTI model. Particularly, in the framework of switched 
LPV controls, it is always desirable to work with a small number of LTI models for control design. This will keep 
the number of controller switching minimal, while avoiding high order dependency on scheduling parameters. In this 
paper, we propose an Adaptive Grid Step Size Determination (AGSSD) algorithm, through which a trade-off between 
the number of local LTI models and the modeling accuracy with guaranteed error bound can be assessed. Finally, each 
LTI model might contain large number of aeroelastic modes; in addition to the rigid body dynamics, which makes 
the control design impractical. A novel model reduction process called Composite Modal Cost Analysis (CMCA) is 
developed to attain a reduced-order model which is better suited for the control design. The essence of CMCA is to 
utilize the collective contribution of each aeroelastic mode throughout the gridded parameter space and retain those 
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modes with significant contribution. 

The approach presented above consists of the following four steps. First, the modal coordinate transformation 
is utilized to convert the original LTI models into the modal canonical forms, where each modal representation is 
decoupled. Second, a modified .#%-norm is used to compute the alignment error between each mode and the modes 
of neighboring LTI models. To better handle the situation where a LTI model has unstable modes, a o-shifted .%4- 
norm, denoted by #(,g_2)-norm, is defined by shifting the imaginary axis of the complex plane to the right by o 
amount, where o > 0 is chosen to be sufficiently large so that the #(,_2)-norm of a transfer function is well defined. 
The (¢_2)-norm of the error transfer function between modes of adjacent LTI models is used as the modal matching 
criteria. This criteria proves to be effective and plays a crucial role in determining the level of alignment, in that 
the smallest value indicates the closest of the two modes and hence should be assigned to the same index order 
throughout the LTI models over the gridded parameter space. Similarly, a o-shifted #.-norm, denoted by Hg_..)- 
norm, is defined as another modal matching criteria. Third, for a given #(g_2)-norm or #(g_..)-norm error bound, 
an AGSSD algorithm is developed to minimize the number of LTI models needed for linear interpolations over the 
gridded parameter space. This algorithm proves to be effective and can significantly reduce the required number of 
LTI models, while providing a balanced trade-off with model accuracy. This result is especially critical for LPV- 
based control design. Finally, a CMCA method is used to obtain a set of reduced-order LTI models for actual flight 
controller design. In summary, the proposed approach not only reduces the total number of LTI models needed for 
covering the entire flight envelop, but also reduces the order of each LTI model so that designing an effective flight 
controller becomes more tractable. 

For ease of presentation, we make use of the NASA Generic Transport Model (GTM) in longitudinal direction 
operating at various flight speeds to illustrate the proposed LPV modeling approach described above. In particular, 
the GTM under consideration consists of a short period mode and 20 aeroelastic modes, and its flight speed ranges 
from Mach 0.5 to Mach 0.88 at a constant cruising altitude. It should be noted that at Mach 0.78, the aeroelastic wing 
is already exhibiting fluttering behavior, since in this condition the first torsional mode becomes unstable. Both time- 
domain and frequency-domain simulations will be used to demonstrate the effectiveness of the proposed modeling 
concept. 

The rest of the paper is organized as follows. Section[2|presents mode alignment procedure that consists of modal 
coordinate transformation and alignment algorithm. A novel AGSSD algorithm is developed in Section[3]to minimize 
the number of local LTI models over the parameter space. A model reduction process using CMCA is presented 
in Section [4] which is used to extract a reduced-order model set containing the dominant aeroelastic modes. Some 
concluding remarks are given in Section|5] 


2. Modal Alignment for Multiple Aeroelastic Models 


This section presents a novel procedure to align aeroelastic modes of multi-input multi-output (MIMO) LTI models 
generated from various operating conditions within the flight envelop. This procedure guarantees the consistency of 
the modal representation for LTI models, which can then be used for LPV controller design. 


2.1. Modal coordinate transformation 


Consider the following collection of LTI dynamical systems, 


i, f X(t) =A'x(t) + Biu(t) _ 
G': { y(t) = Cix(t) i=1,2,---,p. (1) 


where G! represents the aeroelastic LTI model at operating point i, in which x(t) € IR” denotes the state, u(t) € R™ 
the input, and y(t) € R” the output. We assume there are p gridded points in the parameter space within the flight 
envelop, which also corresponds to the number of local LTI models. Furthermore, we assume that each LTT model 
contains m flexible modes. Performing modal coordinate transformation to each local LTI model G! described in 
yields 
i. J X(t) = @'x(t) + B'u(t) _ 
a: { y(t) = x(t) i=1,2,---,p. (2) 
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Figure 1: An illustration of modal alignment process. 


where the system matrices are partitioned as follows, 


; i : : i ph. : Bi . . ; : 
lee 4 i ; aj=| Sh ai oie : , e=[C Se Cn | - (3) 
5 : : 0 J oy : 
Oo -- O Aj, Brn 


Note that Ai denotes a real matrix form for a pair of complex eigenvalues ni = a 25 jbi, where ni denotes the jth 
eigenvalue of the ith LTI model, hence each Ai represents a flexible mode. In this study, we consider two sources 
of modal mismatch when aligning modes from two different aeroelastic LTI models. In general, the sequential order 
of modal frequencies and mode shapes between the two LTI models 4! and Y*, where k € [1,2,---, p| and k Zi. 
That is, the sequential order of block diagonal elements in ./' and .o/* might not be aligned, and this misalignment 
is caused by the notable variations in aeroelastic behaviors from one flight condition to another, such that the order of 
one or more modes is reversed when comparing the two LTI models; see for instance Figure (Ifa). As a consequence, 
performing linear interpolation between Y and Y* will result in an erroneous, mismatched LTI model. Another 
source of modal mismatch is the variations in the plant matrices 4 and @’ that could lead to different normalization 
with respect to each mode [15]. 

Since the aeroelastic dynamics vary continuously over the varying parameters, we adopt the transfer function 
representation to describe the modal dynamics for each given LTI model. 


2.2. Modal alignment algorithm 


From the structure of Y' described in and Qh. the modal transfer function for each mode in ¥! can be expressed 
as 


Gi(s) =C; (sl — Aj) "Bi, J=1,2,---,m; 1=1,2,---,p, (4) 


and the total transfer function for each LTI model ¥ in can be written as a sum of individual modal transfer 
function as follows, 
m m 
G'(s):= L Gils) = LOCA e, f= 12s gp, (5) 
j= j= 
This section is to develop an automated alignment method to match modes with the closest dynamic properties at 
neighboring LTI models, such that all local LTI models will have consistent realizations suitable for LPV control 
design. To this end, we assume that all the local systems are already transformed into the modal transfer function 
form as given in (4p. The modal alignment objective is to identify the closest match of the jth mode at ith grid point 
Gi(s) with the kth modes Gj‘!(s) at the neighboring (i+ 1)th grid point, where j,k = 1,--- ,m. Therefore, once a 
matched mode is identified, the sequential order of modes in Y'*! needs to be permuted so that the jth mode at ith 
grid point is the jth mode at (i+ 1)th grid point. Figure [i{b) shows the aligned aeroelastic modes within each LTT 
model. This framework will be used in Section ]for developing AGSSD algorithm. 
A modified #4-norm is proposed as a means to measure the modal mismatch. Note that the 4-norm for a stable 
jth mode at ith grid point Gi (s) is defined as[27] 


; . ogi 
|| G(s) a= trace(CiPiC; ), 

where Pi > 0 is the unique solution to the following Lyapunov matrix equation, 
AiPi + Pia + BiBi =0. (6) 


Note that the .%-norm is only defined for the stable systems, however, for a flexible wing aircraft, some aeroelastic 
modes can become unstable or flutter at higher Mach number, in which case the #4-norm is not defined. To handle 
this situation, we introduce a notion of relative stability in frequency domain by shifting the imaginary axis to the 
right by o amount, such that in the new shifted complex plane, all modes become stable hence we can compute the 
AG, norm for each modal error transfer function. Mathematically, this frequency shift concept can be described as 
follows, 

Gi(sto)=Ci[(sto)I—Ai]) 'Bi, f= 1,2,---,m; i=1,2,---,p, (7) 


where o > 0 is properly chosen such that it is greater than the largest positive real part of all unstable modes over the 
LTI model set. Therefore, the notion of #(,_2)-norm is defined as follows, 
|| G(s) lo-2=I| G(s + 6) ||2= trace(CiP/C;), (8) 
where Pi > 0 satisfies the following Lyapunov equation, 
(Ai, — o1)Pi + Pi(Ai. — 0)! + Bi Bi =0. (9) 


As mentioned earlier, by shifting the imaginary axis to the right by o, the resulting system matrix Ai — ol will be 
stable so that the HA(g—2)-norm is well-defined. Similarly, we can introduce the #(,_..)-norm of a rational transfer 
function as 

IGj(9)Ilo—-wo= ||Gj(s + loo, (10) 
in which the .#%-norm is defined as 


IGi(s)llo= sup{ Onax(Gi(J@))}, (11) 


where Omax(-) denotes the maximum singular value. 
After transforming all the transfer functions into the shifted frequency domain, the following #¢_2) criterion is 
used to measure the distance between any two modes of the two neighboring grid points, i and i+ 1, 


Mi =||Gi(s) -— Gt" (s)lle—2, i=1,2,---,p—], Jk =1,2,-++,m, (12) 
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Figure 3: GTM with five equally-spaced measurement points along the wing. 


where Mi*! denotes the #(¢_)-norm of the error transfer function between the jth mode at grid point i and the kth 
mode at grid point i+ 1. Note that the smallest value will indicate the two modes have the similar dynamic properties. 
By utilizing for all j and k, we can identify the correct permutation index in Y'*! so that all modes in Y! and Y'*! 
are sequentially aligned. Thus, this alignment process guarantees consistent variations of the flexible modes over the 
entire flight envelope. Algorithm 1| in Fig. 2|illustrates this iterative modal alignment process. 

As mentioned earlier, the proposed approach is applied to the NASA Generic Transport Model (GTM) shown in 
Fig. 3] A set of 18 local LTI models are generated according to various flight speeds ranging from Mach 0.5 to Mach 
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Table 1: GTM models at various Mach numbers. 


ine Local LTI Model G|GP|\G|@P|gG | gi | GP | g@ | gO! gli | gl | gis | gi4 | gl | glo | gl7 | gIB 
ine Mach Number || 0.5 | 0.52 | 0.55 | 0.57 | 0.60 | 0.62 | 0.65 | 0.67 | 0.70 | 0.72 | 0.74 | 0.75 | 0.76 | 0.77 | 0.78 | 0.80 | 0.85 | 0.88 
ine 
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Figure 4: Modal alignment: scrambled (left) vs. aligned (right). 


0.88, as shown in Table[I] We assume that both longitudinal rigid body states and aeroelastic wing states are available 
for measurement. For aeroelastic measurement, as depicted in Fig.|3} we take outputs from 5 equally spaced data 
points for both bending and torsional displacements and their rates. Each local LTI model consists of 20 aeroelastic 
modes and one rigid body short period mode. However, in this study, we focus only on the aeroelastic modes, hence 
we set p = 18 and m = 20. 

To demonstrate the effectiveness of the proposed modal alignment process, the modes within each LTI model 
are intentionally shuffled to generate modal inconsistency among LTI models. Figure [4] illustrates before and after 
applying the modal alignment algorithm, where the left hand side shows the initial scrambled modes and the right 
hand side the aligned modes. It is clear from Fig [4] that, despite the severe modal inconsistency within each LTI 
model, Algorithm | is able to successfully align all the modes sequentially. 

The frequency responses have also been generated to examine modal alignment over the entire speed range. 
Figure [5] shows the frequency responses of mode #13 before and after the execution of Algorithm 1. It is apparent 
that the proposed modal alignment algorithm is able to keep track of all the modes successfully over the entire flight 
profile. It is worth mentioning that .#G_.. criteria have been studied also for mode alignment and produced exactly 
the same result¢2] 


3. Adaptive Grid Step Size Determination (AGSSD) 


With the current advanced computing power and software capabilities, it is not difficult to generate a large number 
of LTI models over the gridded parameter space within the flight envelop, for the purpose of improved model accuracy 


? It is worth mentioning that all flexible modes of the GTM model are complex conjugate. The developed modal alignment algorithm is able 
to align these modes successfully over the entire range of flight speed. However, other scenarios exist where under-damped complex modes could 
become over-damped real modes, and vice versa, as flight condition changes. This problem is called *mode veering phenomena.” The current mode 
alignment algorithm is not designed to handle this general scenario nor to handle the defective systems. However, these will be our future research 
topics. 
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Figure 5: Frequency response of mode#13: scrambled (left) vs. aligned (right). 


for flight control design. However, from the LPV control design point of view, be it switching or non- 
switching [30], a minimum number of local LTI models in the problem setup is always desirable, because this would 
minimize the number of local controllers needed for control design, reduce the number of Linear Matrix Inequalities 
(LMIs) required for controller synthesis, and eliminate numerical issues during the control design. Now, the challenge 
is how to reduce the number of LTI models without sacrificing model accuracy. In this section, we propose a novel 
adaptive algorithm to adjust the grid step size with prescribed error bound. Again, we utilize the notions of H,_2)- 
and #(_..)-norm defined earlier as metrics to measure the modeling error. 

To this end, we assume all the LTI models in the gridded parameter space have been aligned by following the 
modal alignment algorithm, Algorithm 1. Let || G'*(s +.) ||2 denote the Hg—2)-norm of the modeling error between 
G' and Y'**, ie. . . . 

IF (s) o-2=IF'() 905) Ilo-2 - (13) 


Let 6 > 0 be a prescribed HAG g—2)-norm error bound, then all the LTI models between ¥ ‘and G'+*, k > 1, can be 
eliminated, if 


|F'*(s) |lo-2< 6 and ||G**!(s) ||o_2> 6. (14) 


In other words, if the conditions in are satisfied for a given error bound 6, all the LTI models between indexes i 
and i+k can be considered redundant and hence eliminated without much impact on the overall model accuracy. The 
detailed AGSSD process is described in Algorithm 2. Similarly, we can also utilize the #(,_..)-norm to compute the 
modeling error, i.e. 


WF" (s) o-0=IF'(s) —F'*(s) [oe (15) 


and the same AGSSD algorithm can also be applied to attain the reduced number of LTI models for a given 6 error 
bound. Note that the choice of 6 is critical, for it dictates the number of retained LTI models at the end of AGSSD 
algorithm. To illustrate this, both #_2)- and #(,_..)-norm are used in AGSSD algorithm to process the 18 LTI 
aeroelastic GTM models given in Table 1. Figure |[6|shows the number of retained LTI models as functions of H,_2) 
and Hg —c) error bound. As expected, these two norms render different results. However, both norms show that as 
the requirement on the error bound becomes tighter, additional LTI models will need to be retained in order to meet 
the prescribed level of model accuracy. For instance, for #(¢_2)-norm error bound of 300, we only need to retain 
eight LTI models out of the original eighteen models to cover the entire flight regime. However, it is still unclear as 
to what that specific 6 error bound means in assessing the level of model mismatch between any two LTI models. 
To address this, for a given error bound, we need to examine both the time domain and frequency domain responses 
from the two LTI models and correlate the level of model mismatch to the error bound. Through this exercise, we can 
determine the acceptable 6 error bound and hence the number of retained LTI models in the gridded parameter space 
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algorithm 2: Adaptive grid step size determination (AGSSD) algorithm. 


e Set H@g_)-norm (or #G~.-norm) error bound 6. 


e Leti< 1 and Ind < [1]; and given 6 > 0 and number of LTI models p. 


while i < p—1do 

e Letk <1. 
while k < p—1do 

Calculate || Y'*(s) ||2_, using (or || Y'*(s) ||2_... using (15)). 
if 5 >||9'*(s) |[G_o (or 6 >|] F"*(s) ||. ) then 

kek+l1 


cy 


else if k = 1 then 


Ind < [Ind i+k], ic-it+k; 
else 
Ind < [Ind i+k—l], ic it+tk—-1, ke1; 
end 
if i= p—1 then 
Ind < [Ind i+k); 
Break 
else if i+-k > p then 
Ind < [Ind i+k-—1]; 
Break 
end 
end 
ifi+k > p then 
Break 
end 


end 


to be used for subsequent model reduction. Figures [7|to[10]show the time domain and frequency domain comparisons 
of the linearly interpolated LTI models, © (Mach 0.62) and Y!3 (Mach 0.76), at two levels of error bound, 5 = 300 
and 6 = 600, based on HA g—2)-norm. As shown, when 6 = 300 the responses from interpolated models match very 
well with those from the exact model, whereas when 6 = 600 the responses of the interpolated model shows a large 
deviation from the actual model. 

On the other hand, when utilizing #_.-norm as a metric, at 6 = 300 we see that we only need 5 LTI models to 
cover the entire flight envelop. However, for validation purpose we perform a linear model interpolation between 4! 
(Mach 0.5) and ¥!? (Mach 0.75) to attain the interpolated model Y 7 (Mach 0.65). Figure|11]shows the time-domain 
simulations of the interpolated Y’ model versus the actual Y’ model, and it shows a large deviation. This indicates 
that .#G_. error bound of 300 is not acceptable, though it renders fewer number of LTI models. Furthermore, if 
we choose .#G_. error bound of 200, we will need 6 LTI models, saving one LTI model compared to the case with 
Ag —2 error bound of 300. Figure|12|shows the validation of the interpolated Y> (Mach 0.6) derived by interpolating 
between Y! (Mach 0.5) and Y* (Mach 0.67). Similarly, Fig. [13]shows the validation of the interpolated Y!° (Mach 
0.72) by interpolating Y* and Y'? (Mach 0.76). These time-domain simulations reveal that .7_.. error bound of 200 
renders an acceptable result. 
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Figure 6: Number of retained LTI models as function of error bound. *: #(g_2) norm; 0: H#(g—..) norm. 
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Figure 7: Time-domain simulations as function of error bound (Y 6 Mach 0.62). 
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Figure 10: Frequency-domain responses as function of error bound (¥!3, Mach 0.76). 
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12 


LTI#5 Interpolated via LTI#1 & LTI#8 
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Figure 12: Time-domain simulations of the interpolated Y° model vs. actual model at .%_.. error bound of 200. 
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Figure 13: Time-domain simulations of the interpolated Y!° model vs. actual model at .#_.. error bound of 200. 
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4. Model Reduction 


Because of the practical limitations on control actuation bandwidth and the needs to lower computational burden 
in actual implementation, model reduction is often an essential part of modeling and control effort for any physical 
systems of large dimensions. Following the two algorithms presented in the previous sections, we have attained a 
reduced number of LTI models with all their modes properly aligned within each LTI model. However, each LTI 
model still contains very high number of flexible modes. Hence, in this section we introduce the notion of Composite 
Modal Cost Analysis (CMCA), through which we can attain a reduced-order model that is better suited for flight 
control design. 


4.1. Composite modal cost analysis 

The basic idea behind Modal Cost Analysis (MCA) is to examine the contribution of each flexible mode 

to the mission objectives in a control system. A metric of modal contribution can be calculated in terms of output 

covariance, from which contribution of each mode is evaluated and ranked from high to low. This approach is used 

to derive a reduced-order model from a full-order model by neglecting less significant modes. Recall the modal 
coordinate representation of a LTI model ¥! described in (2), subjected to the disturbance input, 

gi - { x(t) = @'x(t) + Blu(t) + F'w(t) 


yh =ex0) i=1,2,---,p, (16) 


where w(t) € R”” is the random disturbance input with intensity W > 0 and the system matrices (./', Z',@') are 
block matrices as given in (). Similarly, the matrix Y' can be partitioned accordingly as follows, 


Dh 
If o/' is Hurwitz, then the open-loop output covariance cost for 4! is given by 
Yi:= trace(C'P'S"), i=1,2,---,p, (17) 
where P’ > 0 is the controllability Gramian matrix satisfying the following Lyapunov equation, 
Pid’ to'P+9WF =0, i=1,2,---,p. (18) 


Let vi denote the jth modal contribution of ith LTI model 4! to the output covariance, then noting the modal block 
partitioning of the matrices (.e/', ',@'), we can deduce that 


m 
Yi=)ivi; vi:=trace(C\P'C}),  i=1,2,---,p, (19) 
j=l 


where Cr ! is a 2 x 2 output covariance matrix corresponding to the jth mode in Y’. This indicates that the total 


output covariance cost for Y! can be expressed as a collection of its individual modal contribution. Therefore, we can 
compute the modal cost of each mode v; from and rank its contribution from high to low as 


[vil > v5] > ++ > [Vil (20) 


where vi is the most critical mode and v’, the least critical mode among vi, j=1,---,m. Note that the modal cost vi 


can be a negative value, which indicates that this particular mode is in fact helping to reduce the total cost Y', however 
the total cost is non-negative. Since MCA involves solving the Lyapunov equation (18), it can only be applied to 
the stable aeroelastic modes for model reduction. Hence, for conventional MCA to work, one must first identify 
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Table 2: Composite modal cost of the stable aeroelastic modes. 


ine Mode index 
ine Total modal cost 


jas Vania] tena ae al Gere | ate | 8) See te leo eee | ee |e a | 
ine 


468.36 | 357.93 | 314.84 | 161.72 | 133.46 | 99.66 | 82.56 | 53.19 | 33.37 | 28.87 | 23.10 | 13.13 | 11.11 | 7.28 | 6.00 | 3.49 | 1.36 | 0.55 


Modal cost 
nN 
8 


LTI models 


Figure 14: Modal cost of the stable aeroelastic modes. 


and decompose the modes into stable and unstable modes, and perform MCA only to the stable modes for model 
reduction. It should be noted that the unstable aeroelastic modes and the aircraft rigid body modes are to be retained 
by default in the reduced-order model. 

To perform the MCA for multiple LTI models covering a wide range of flight regime, we utilize CMCA, in which 
the collective contribution of each mode is summed up throughout gridded parameter space and its contribution ranked 
from high to low. In other words, if 7%; denotes the collective contribution of jth mode over all grid points, then its 
composite modal cost is given by 


P o 2 ae Pee A 
M; = ¥ vi; vi :=trace(CiP'C;), j=1,2,-+-,m. (21) 
i=! 


Similarly, we can compute the composite modal cost for each mode throughout the gridded parameter space and rank 
its contribution as 


|M@| >= || >--->|Gnl, (22) 


where .@ is the most contributing mode and .%,,, is the least contributing mode among .@;, j = 1,2,--- ,m. Next, 
we can utilize the proposed CMCA to each of the 20 aeroelastic modes throughout the 18 LTI models, and the results 
are given in Table[2] Since there are two unstable modes (1st bending and Ist torsion) at higher flight speeds, they are 
omitted from CMCA computation, hence only 18 modal costs are shown in Table [2| Figure[14|shows the modal cost 
of Vj for i,j = 1,--- , 18, from which the composite modal cost is derived. 

To better handle the cases with unstable modes, we extend the notion of CMCA by introducing a o-shift transfor- 
mation to all the modes throughout Y', i = 1,--- ,18, so as to ’stabilize” all the LTI models in o-shifted coordinates. 
Table 3] lists the composite modal cost for all 20 aeroelastic modes after o-shift, whereas Figure [15] shows the in- 
dividual o-shifted modal cost. Careful examination of Tables [2] and 3] reveals that the most contributing modes are 
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Table 3: Composite modal cost of all aeroelastic modes with o-shift. 


ine Mode index 
ine Total modal cost 
ine 


2 4 6 3 8 I 7 10 5 18 13 9 11 17 16 12 15 14 20 19 
322.84 | 299.39 | 188.44 | 177.54 | 158.35 | 136.86 | 105.93 | 104.78 | 67.30 | 62.72 | 38.91 | 34.45 | 27.06 | 24.23 | 17.73 | 12.26 | 8.97 | 7.56 | 3.20 | 1.48 


Modal cost 


LTI models 


Figure 15: Modal costs of all aeroelastic modes after o-shift. 


consistent and that the summation of the first four modal costs from Table 2] or first six modal costs from Table 
amounts to more than 75 % of the total modal cost. Note that the two additional modes italicized in Table [3]are pre- 
cisely those omitted unstable modes from Table [2| Therefore, we can choose six aeroelastic modes out of 20 modes 
throughout Y', i = 1,--- ,18, to form a reduced-order model suited for control design. Table|4]shows an example of a 
reduced-order model containing 6 most significant modes at Mach 0.88. 


4.2. Model validation 


To validate the proposed model reduction process, we present the root-locus, and the time- and frequency-domain 
comparisons between the full-order model and the reduced-order model at various flight conditions. For this study, 
the reduced-order model is consisted of one rigid body short period mode and six aeroelastic modes chosen through 
the CMCA procedure presented earlier. 

Figure[16|shows the root-locus of the full-order models and reduced-order models over the entire flight envelope. 
This figure demonstrates the smooth transition of the poles over all grid points, which also verifies alignment of 
the reduced-order models. It is clear that the unstable dynamics are kept intact in the reduced-order models and, as 
expected, the most contributing modes are close to the origin. 

Furthermore, a series of time-domain simulations for full- and reduced-order models at Mach 0.7 (4?) and Mach 
0.8 (Y'®) are conducted, and they are shown in Figs. [17] and respectively. Except at the transient region (for 
torsional displacement and torsional displacement rate) where high frequency contents dominate, overall the reduced- 
order model is able to successfully capture both rigid body dynamics and wing-tip aeroelastic behaviors of the full- 
order model. Similarly, Figure [19] and [20] show the frequency response comparisons between the full- and reduced- 
order model at the two flight conditions. The responses are taken from the disturbance input to the the wing-tip 
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Table 4: A reduced-order model with six aeroelastic modes at Mach = 0.88. 


ine Mode ID Frequency (rad/sec) | Damping ratio 
ine Ist bending (1) 10.5202 -0.0926 
ine 2nd bending (2) 11.27 0.3158 
ine 3rd bending (3) 17.0996 0.3451 
ine Ist torsion (4) 18.3889 -0.1268 
ine 2nd torsion (6) 35.4116 0.0291 
ine 4th bending (8) 68.4175 0.0251 
ine 
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Figure 16: Root-locus across flight envelope. 
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Figure 17: Time-domain simulations at wing tip, Mach = 0.7. 


bending displacement and pitch rate, respectively. These figures show that the reduced-order model is capable of 
capturing the dynamics of the full-order model at the frequency range of interest. 


5. Conclusions 


In this paper, we developed a novel modal alignment algorithm to sequentially order the flexible modes for a 
collection of large dimensional LTI models sampled over the gridded parameter space within the flight envelop. The 
modal mismatch criteria based on the o-shifted .#%- and .#%2.-norm were defined and utilized to align modes of similar 
dynamic behaviors at neighboring LTI models. An Adaptive Grid Step Size Determination algorithm was developed 
to minimize the number of local LTI models needed to cover the entire gridded parameter space with guaranteed error 
bound. This step is especially critical when designing LPV-based flight controllers, in which a problem setup with 
a small number of local LTI models is highly desirable. Finally, we proposed the concept of composite modal cost 
analysis and utilized it to attain a reduced-order model that captures essence of the full-order model. Throughout 
this paper, we used the NASA GTM aeroelastic aircraft models to illustrate the developed algorithms and the model 
reduction process. The analysis and simulation results demonstrated the effectiveness of the proposed concept. 
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Figure 18: Time-domain simulations at wing tip, Mach = 0.8. 
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